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Quark-Model Baryon-Baryon Interaction Applied to the 
Neutron-Deuteron Scattering (III) 

Breakup Differential Cross Sections 

Yoshikazu Fujiwara and Kenji Fukukawa 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

The low-energy breakup differential cross sections of the neutron-deuteron (nd) scattering 
are studied by employing the energy-independent version of the quark-model baryon-baryon 
interaction fss2. This interaction reproduces almost all the breakup differential cross sections 
predicted by the meson-exchange potentials for the neutron incident energies E n < 65 MeV. 
The space star anomaly of 13 MeV nd scattering is not improved even in our model. Some 
overestimation of the breakup differential cross sections at E n = 22.7 - 65 MeV implies that 
' systematic studies of various breakup configurations are necessary both experimentally and 

' theoretically. 
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§1. Introduction 



The three-nucleon (3iV) system is a good place to study the underlying nucleon- 
nucleon (NN) interaction, since many techniques to solve the system exactly are 
well developed nowadays Ample experimental data are already accumulated 
especially for the low-energy neutron-deuteron (nd) and proton-deuteron (pd) scat- 
tering and extensive studies to detect the 3iV force have been carried out based on 
CN ! the modern meson-exchange potentials,®® and more recently, on the chiral effective 

field theory.® ® Most of the researches to such a direction are concerned with higher 
£f~^ ■ energies than 100 MeV for the nucleon incident energy E n in the laboratory (lab) 

system, since the 3N force effect is expected to be revealed more prominently than 
in the low energies.® On the other hand, the discrepancies of various 3iV observables 
between the theory and experiment in the E n < 65 MeV region, are not resolved 
even by the recent accurate treatment of the Coulomb forced - ® This is particularly 
true for the nucleon-induced deuteron breakup processes. It is therefore worth while 
reexamining the NN interaction itself if the present-day realistic force is the most 
appropriate one to start with. 

In previous papers ES'EJ referred to as I and II hereafter, we have applied the 
quark-model (QM) baryon-baryon interaction fss2 to the neutron-deuteron (nd) elas- 
tic scattering. This interaction model, fss2, describes available NN data in a com- 
parable accuracy with the modern meson-exchange potentials.^ By eliminating the 
inherent energy dependence of the resonating-group kernel, fss2 was found to yield 
a nearly correct triton binding energy, the S-wave nd scattering length, and the 
low-energy eigenphase shifts without reinforcing it with the three-body forceP®~E® 
The predicted elastic differential cross sections have sufficiently large cross section 
minima at E n = 35-65 MeV and 8 cm = 130° - 135°, in contrast to the predictions by 
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the standard meson-exchange potentials.^ 1 The so-called A y puzzle at low-energies 
E n < 25 MeV is largely improved in this modelP^ In this paper, we continue these 
studies by examining the 3iV breakup processes with various decaying kinematics for 
the energy range E n < 65 MeV. The main motivation is to find if the quite different 
off-shell properties, originating from the strong nonlocality of the QM baryon-baryon 
interaction, give some influence to the 3iV breakup differential cross sections. In 
contrast to the elastic scattering amplitude, the breakup amplitude covers a wide 
momentum region of the three-body phase space. It will be found unfortunately 
that the fss2 gives predictions similar to the meson-exchange potentials and does 
not improve much the discrepancies between the theory and the experiment. 

The organization of this paper is as follows. In §2.1, the formulation of the 
breakup differential cross sections is given in terms of the direct breakup ampli- 
tude. Various kinematical configurations for the three-body decay are introduced 
in § 2.2. A minimal description of the three- nucleon breakup kinematics is given 
in Appendix A. The isospin factors for the breakup amplitudes are derived in Ap- 
pendix B. The comparison with the experimental data is presented in § 3 for energies 
E n = 8, 10.3, 10.5, 13, 16, 19, 22.7 and 65 MeV. The difference from the predic- 
tions by meson-exchange potentials are discussed in detail. We close this paper with 
a summary of this series of investigations in § 4. 

§2. Formulation 

2.1. Breakup differential cross sections 

Following the notation of Refs.[Tj), I and II, the three-body breakup amplitude 
is given by 

U \(f>) = (1 + P)T\<f>) = (1 + P)tQ\4>) . (2-1) 

In order to derive the breakup differential cross sections, we start from the Fermi's 
golden rule 

2tt 1'°° 
dN = —\(pq\U \4>)\ 2 p 2 d P 5(E-E pq ) q 2 dq dp dq , (2-2) 

and divide it by the incident flux j = (3% /2M)/(2vr) 3 . Here, E pg = (h 2 /M)(p 2 + 
3q 2 /4), M is the nucleon mass, and qo is the incident momentum related to the 
energy, E = (3h 2 /4M)qo 2 + Ed, in the center-of-mass (cm) system. We obtain in the 
cm system 

jdpdqdq 3h 2 q J 

(27r)4 (^)lv ! 5: Kwr|c/ ^ )o12 ' (2 - 3) 



d 5 a 
dp dq dq 
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where Eq. 1(2.880 is used to perform the p-integral. In Eq. (|2-3p . r = r a r T is 
the spin-isospin quantum numbers in the LS-couphng scheme and the subscript 
in the matrix element implies the on-shell condition |p| = po = \J (3/4)(gjvf 2 — Q 2 ) 
with qM = v qo 2 — K-d 2 - Here, is related to the deuteron binding energy \e^\ 
through | | = (3h 2 / AM)Kd 2 ■ In this paper, we use the notation r a = {s\)SS z and 
r r = (t\)\T z to specify the quantum numbers in the -LS-coupling scheme: i.e., 

|p,<?;123) =J3b,g,7> <7|P,9;123) , 

7 

(p,q;123| 7 > = [Y {xe)L (p,q) £ r<r (12,3)] ^ r,r T (12, 3) , (2-4) 

with 7 = [(A^L-To-] JJ«; r T , and £r CT and ??r T being the three-particle spin and isospin 
wave functions, respectively. In Eq. (|2-4p . I^a^lmCP) <?) = [^A(p)^K<z)]LAf are the 
angular functions. For the initial state, we use channel-spin representation as for 
the elastic scattering. We take the sum of Eq. (|2-3|) over all the spin and isospin 
quantum numbers and divide by the initial spin multiplicity 6. The selection of the 
detected particles in the final state is controlled by the isospin projection operator 
O r , the explicit form of which will be specified later. The breakup differential cross 
sections of the nd scattering are therefore calculated from 

r ScScz 

(2-5) 

Let us first consider the spin-isospin sum I = J2r KPQ-TKl+P)/)^ by neglecting 
the initial spin quantum numbers for the time being. The effect of the permutation 
P (123) Q in (1 + P) = J2 3 a=1 P {l23) a is defined by 

(pq\P(123) a f) = P(123) a f(P, 0) = f(Pa, la) , (2-6) 

if the function f(p,q) does not contain the spin-isospin degree of freedom. In fact, 
we should use 

(pqr\P {123 ) a f) = (r\P$l"f{ Pco q a )) , (2-7) 

where Pq%L is the permutation operator in the spin-isospin space and the bra-ket 
notation is used for the spin-isospin degree of freedom. Using these notations and 
the completeness relationship in the spin-isospin space, ^2 r \ r)(r\ = 1, we find 

I = E {f(Pa,Qa)\P^)^P§S) P \f^Qp)) ■ (2-8) 

Here we separate the a, (3 sum into the diagonal part (a = f3) and the off-diagonal 
part In the off-diagonal part, we specify a and (3 by the cyclic permutations 



*' In the following, we cite equations of the previous paper I (or II), with adding I (or II) in 
front of the equation number. 
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Of = {P[x2*T a O t (P[ 123) Y (= O^) , (2-11) 



of (123) ((a/?7)=(123)-cyclic). For these terms, the a-/3 term and the /3-a term are 
complex conjugate to each other. Thus we obtain 

3 / 

J = £</(p«,OI/(p«,9a)>+2 E ^(/(Pa.gJI^S)!/^.^)) , (2-9) 

a=l (a/97) 

where ^ implies the sum over the three cyclic permutations of (a/37) = (123). 

The extension to I = Ylr \(pQP\Ot(^ + P)f)\ 2 i incorporating the isospin pro- 
jection operator O r , is rather easy. Here, O r is specified as 

opp = 1 + t z (1)1+t z (2) Qnn = 1-^(1)1-^(2) 
2 2' 22' 
= 1 + t z (1)1-t z (2) = 1-t z (1)1 + t z (2) 

2 2 2 2 ' 1 J 

depending on the species of particles 1 and 2 detected. We use O r 2 = O t and the 
notation (pqr\f) = fripiQ)- Then, by defining 

3-a / \/9 
(123) ) O t [P(\23) J 

we obtain 

3 

1 = E E /? (p«> ««) <^icrm /r(p„, o 

+2 E ERe{/i(p a ,qJ(f|P ( CT 123) 0^|r)/ r (p / 3,q / 3)} . (2-12) 

(a/37) f ,r 

The spin-isospin factors in Eq. (|2-12p are calculated by separating the spin-isospin 
state \r) to the spin and isospin parts, \T) = \r a )\r r ). We find 

(r\o?*\r) = 5§ tS 6j iS (r T \o? a \r T ) , 

(r \p^ m of\r) = s St8 (-i) 1+s xl (r T \of\r T ) , (2-13) 

where r a = (s±)SS z and F T = (t\)\T z , and Eq. (jB-ip is used for the spin part. We 
also extend the definition in Eq. (|B- 1|) for the spin part to the isospin part as in Eq. 
(IFT21) . Using the definition of Xl (a/3) in Eq. (IF2l) . we can write the matrix elements 
in Eq. (I2d3j) as 

(r\or\r) = 8 §tS 5j >s xi< t aa) , 

(r \P[ m of\r) = s StS (-i) 1+s xl (-i)^xz^ , (2-w) 

for (a/?7) = a cyclic permutation of (123). Thus we find 

3 

1 = E E ^ f r ( p °" q ^ 

a=l p p 
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+ E Y. 5 ~s,s(- 2 ) X lsX~^ ] Re{r r ( Pa ,q a )fr(- Pf} ,q p )} • (2-15) 

(a/37) r,r 

where the generalized Pauli principle (— l) s +'+ A = — 1 is used for the two-nucleon 
part of fr{—pp,qp)- The isospin factors A-j ^ are explicitly given in Appendix B. 
We assign the direct breakup amplitude to fr{p,q) in Eq. (|2-15p through 

fr,s cScz (P,Q) = -(2vr) 2 (pqr\T\ct> qo -S c S cz )o , (2-16) 

with T = tQ. The partial wave decomposition is given by 

7,1', J Jz M 

x Y,(t'™'S c S cz \JJ z ) Y e 1 m ,(q ) , (2-17) 
m' 

where the prime on the sum implies that we take all the orbital angular momen- 
tum sum for the LS coupling scheme of 7; i.e., the sum over only {\t)L with 
7 = [(X£)Lr a ] JJ Z ; r r . (Note the extra (47r) factor for the scattering amplitude.) 
It is convenient to define 

fi(lS B )J _ sr^ fi(t'S' c )J fJ ,„ 1R v 

~~ 2^ ^/n J(l'S' c )(tSc) ' \ Z ' L °) 

by the solutions, Qi^ = Pi 2 Wiq^yfcJ^(pi,q^,j\Q\'ip}, of the basic AGS equation 
in Eq. 1(2.59), and the elastic scattering amplitude f^/ S i^£ S \ in Eq. 1(2.92). The 

partial- wave amplitude for the direct breakup, f!^,g )(?)E3 is expressed as 

= £<poim*w/j»oir) sm-^=q%7 c)j > ( 2 - 19 ) 



if we use the the spline interpolation for a particular value of q. For the practical 
calculations, it is convenient to adopt a particular coordinate system with q = e z 
in Eq. (|2-17|t . Then the basic direct breakup amplitude in the spin-isospin space is 
calculated from 

fr,s c s cz (P,Q) = v 7 ^ Yl 4$s c ) ~S Z )SS Z \JS CZ ) 

7,1', J 

x?(£'0S c S cz \JS cz )Y {Xi)L{Scz _ Sz) ( P ,q) , (2-20) 



*-* This corresponds to the direct term of f^fg/g in Eq. 1(2.92). 
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with g = y/2t' + 1. The differential cross sections in the cm system are given by 



d 5 a _3p q 2 l 
dp dq dq 4 q$ 6 



- y 

Sc-fScZ 



E E <%,S ^,sX~f a) f*^ ScS JPa, q a )fr,S cScz (P a , Qc 



a=l 



r,r 



+ E E^(- 2 )*l^ t M) Re{/i 5c5cz (p Q , g J/ rA5cz (-p„^)} 

(a/37) r,r 



(2-21) 



The breakup differential cross sections in the lab system, (d 5 a / dk\dk,2d S), are 
specified by the two directions k\, k2, and the energy S measured along the locus of 
the E1-E2 energy plane. They are obtained from Eq. (|2-2ip by a simple change of 
the phase space factor 1 ' 



3W 



4 qo 



Plab 



M3kik 2 
V 2q 



k\ 

, COS 6 2 + — COS 012 

k 2 k 2 



— COS U\ + — COS 012 

k\ k\ 



-1/2 



(2-22) 



The details of the three-body kinematics are summarized in Appendix A. 

2.2. Three-nucleon breakup kinematics 

Assuming that we detect two outgoing particles 1 and 2, the breakup differential 
cross sections are specified by two polar angles 0\, 82, and a difference of azimuthal 
angles </>i2 = <f>i — $2, in addition to the energy S determined from the kinematical 
curve (S"-curve) in the E\ and E% plane. The starting value of the arc length S = is 
quite arbitrary and we follow the convention by the experimental setup. In Appendix 
A, we have parametrized the locus in the k\-k2 plane with an angle 9, and the starting 
point S = is uniquely determined by specifying 6 st . We also assume that the beam 
direction of the incoming particle is the z axis and set = 7r,® which determines 
the x-axis. 

It is customary to classify the three-body breakup kinematics into the following 
six categories based on the classical (or geometrical) argument 



1. The quasi- free scattering (QFS): one of the nucleons in the final state is at rest 
in the lab system. 

2. The final-state interaction (FSI): the relative momentum of the two outgoing 
nucleons is equal to zero. 

3. The collinear configuration (COLL): one of the outgoing nucleons is at rest in 
the cm system, and the other two have momenta back to back. 

4. The symmetric space star configuration (SST): the three nucleons emerge from 
the reaction point in the cm system, keeping equal momenta with 120° relative 
to each other and perpendicular to the beam direction (on the x-y plane in the 
cm system). 
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5. The coplanar star configuration (CST): the same with the symmetric space star 
configuration, but with the three momenta lying on the reaction plane. 

6. The non-standard configuration (NS): the other non-specific configurations. 

These are mathematically distinguished by particular values of the lab momentum 
k a , the relative momenta p a and q a , etc., and provide a rough guidance to which 
portion of the two-nucleon i-matrix is responsible at the final stage of the reaction, 
according to the structure of the direct breakup amplitudes in Eq. (|2-19|) . For exam- 
ple, Ref. [1]) argues that the first Born term of the QFS with k a = is approximately 
a product of an on-shell two-nucleon i-matrix and the deuteron state at zero mo- 
mentum. It is known that the 3iV force effect is rather small for the QF condition. 4 ' 
On the other hand, the collinear configurations with q a = are expected to be sen- 
sitive to the 3N force intuitively, and the experimental study by Correll et alf^ was 
carried out to study the effect of the 3N force intensively in the reaction H(ci, 2p)n 
around these configurations at the deuteron incident energy Ed = 16 MeV. Further- 
more, the FSI is characterized by p Q = 0, for which the half off-shell t-matrix in 
Eq. (|2-19p generates a large peak corresponding to the 1 5'o positive-energy bound 
state near the zero-energy threshold. Although the height of the peak is influenced 
by the background amplitude Qi^, the FSI peak is usually nicely reproduced. The 
disagreement with the nd data is reported at the early stage for the SST configura- 
tion, which is still an unsolved problem called space star anomaly!^ 1 It should be 
noted, however, that the disagreement between the theory and experiment is also 
seen in some other coplanar star and non-standard configurations, for which off-shell 
properties of the two-nucleon i-matrix is expected to play a role in different ways. 
We will examine these case by case in the next section. 

§3. Results and discussion 

3.1. H(d,2p)n reaction at E^ = 16 MeV 

It is important to take enough number of discretization points and partial waves 
to get well converged results, especially for the breakup differential cross sections. 
In this paper, we take ni-n2-n3=6-6-5 in the notation introduced in § 3.1 of I, unless 
otherwise specified. This means that the three intervals, [0, </jw/2], [<?m/2, (v / 3/2)q , m] 
and [(v / 3/2)<?m ) Qm], are discretized by the six-point Gauss Legendre quadrature for 
each, and the total number of the discretization points for q is 35 (= 6x3 + 6x2 + 5). 
The three-body model space is truncated by the two-nucleon angular momentum 
J max , which depends on the incident energy of the neutron. We find that I max = 3 
is large enough for E n < 19 MeV. The Coulomb force is entirely neglected in the 
present calculations. 

We first investigate the Correll et al.'s experiment^ i.e., H(d, 2p)n reaction 
with the deuteron incident energy E^ = 16 MeV. This corresponds to the nucleon- 
induced breakup reaction of the nucleon incident energy of 8 MeV. We generate 
the direct breakup amplitude for the deuteron incident reaction by adding an extra 
phase factor (— 1) to each term of Eq. (|2-2U|) . corresponding to the change from q 
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Fig. 1. Breakup differential cross sections for the reaction d(n,2n)p with E n — 8 MeV. The 
experimental data are the deuteron incident reaction, H(d,2p)n. with the energy Ed = 16 
MeV. 113 



to — q in Eq. (|2-17[) 5] The decay kinematics for the deuteron incident reaction is 
discussed in Appendix A. The breakup differential cross sections for the d(n, 2n)p 
reaction with E n = 8 MeV are compared with the Correll et al.'s data in Fig.[H 
with respect to two collinear (COLLI, COLL2) and two non-standard (NS1, NS2) 
configurations. The starting point S = is chosen as the collinear points or the 
nearest point, as is discussed in Appendix A. The dashed curve, the solid curve, and 
the bold solid curve correspond to the S+D (i.e., 3 5i + 3 -Di and x 5o only), / max = 2, 
and I max = 3 cases, respectively. The solid curves almost overlap with bold curves 
and I max = 2 is actually good enough at this energy. We also see that even the 
restriction to the S+D model space is not too bad. If we compare our results with 
meson-exchange predictions in Ref. [T]), we find that they are very similar to each 
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other. The calculated values are somewhat too small especially in COLLI, COLL2 
and NS1, although to less extent for the meson-exchange predictions. There is no 
exact collinear point in the case of NS2, and the best agreement with experiment is 
obtained in this case. 

3.2. d(n,2n)p reaction at E n = 10.3 MeV 

The breakup differential cross sections for the d(n, 2n)p reaction at E n = 10.3 
MeV are displayed in Figs. [2] and [3l together with the experimental data by Gebhardt 
et alP^ 1 The figure number, fig. 5 etc., in each panel corresponds to the original 
one in Ref. ll9|) . The large two peaks seen in fig. 5 - fig. 13 are the np final state 
interaction peaks with p l ~ on the lower S side and those with p 2 ~ on the 
higher S side. On the whole, the comparison with experiment gives fair agreement, 
but some discrepancies found in Ref. ll9p still persist. In Ref.ll9p. the experimental 
data are compared with the solutions of the AGS equations in the W- method, using 
a charge-dependent modification of the Paris potential. Their results and ours are 
strikingly similar to each other, sharing the same problems for the detailed fit to the 
experiment. The peak heights for the p 2 ~ final state interaction peaks are not 
precisely reproduced in fig. 5, fig. 7, fig. 8 and fig. 10, probably because we did not 
take into account the charge dependence of the two-nucleon interaction. In fig. 9, 
our result is worse than the theoretical calculation in Ref. 19). The collinear point 
is realized at S = 4 MeV in fig. 11, at S = 7.5 MeV in fig. 13, at S = 6 MeV 
in fig. 14, and at S = 5.6 MeV in fig. 15. In fig. 11, we have obtained a smooth 
curve around the collinear point, just as the theoretical calculation in Ref. ll9p . The 
breakup differential cross sections around the collinear points are well reproduced. 
In fig. 16, the flat structure between S = 3 - 7 MeV is just the same as the theoretical 
calculation in Ref. 19). The experimental data of Ref. ll9p for the symmetric space 
star configuration is plotted in the first panel of Fig.[6l Here again, we have obtained 
very similar result with the theoretical calculation in Ref.ll9p. 

3.3. Quasi-free scattering 

We show in Fig.[4]the breakup differential cross sections for the quasi-free scatter- 
ing at energies £i a b = 10.5 - 65 MeV. We find some deviation from the experimental 
data at the peak position for all the energies. Detailed investigation of the Coulomb 
effect in Ref.[8]) has revealed that this overestimation at the peak position is reduced 
to some extent. However, the reduction is probably not large enough except for 
the case of E^b = 13 MeV. Figure 8 of Ref.[5]) implies that this reduction is energy 
dependent. The large overestimation at -Ei a b = 19 MeV may not be resolved only by 
the Coulomb effect. The direct incorporation of the Coulomb force is necessary for 
our QM interaction. In Fig.JU we can see that the roles of higher partial waves are 
important for higher energies. The partial waves up to I max = 3 is clearly necessary 
for £"i a b = 19 MeV. For the energies, E^ = 22.7 and 65 MeV, we need more partial 
waves up to I max = 4. 
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Fig. 2. Breakup differential cross sections for the reaction d(n, 2n)p with E n = 10.3 MeV. The 
experimental data are taken from Ref.ll9[) with the same figure number, fig. 5 etc., in each panel. 
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3.4. Final state interaction 

Four examples of the breakup differential cross sections for the np final state 
interaction are displayed in Fig.[5j Almost all the data are for the d(p,pp)n reaction. 
In the 10.5 MeV and 13 MeV cases, the pd data are shown with open circles, while 
the nd data with bars, some dots and diamonds. The lower peaks are the np final 
state interaction peaks with p l ~ 0, while the upper with p 2 ~ 0. Here we find that 
the higher peaks are slightly too small. The Coulomb correction increases the peak 
height a little,^ and improves the fit to the experiment to some extent. We probably 
need more careful treatment of the charge dependence of the NN interaction, just 
as in the previous 10.3 MeV case. We also see that the minimum point at S = 11 - 
12 MeV for the = 16 MeV reaction is too low. The three-nucleon force might 
be necessary to increase the cross sections and get a good fit to the experiment.^ 

3.5. Symmetric space star configurations 

As is well known, a large discrepancy appears in the breakup differential cross 
sections in the symmetric space star configurations. This is seen in Fig.(6j where 
our results in the various model spaces are compared with the nd and pd data. 
A strange thing is that any theoretical calculations of the d(n, nn)p reaction at 
-E-iab = 13 MeV deviate largely from the old and new nd data,^'^'^ although the 
deviation is not much for the 10.3 MeV, 19 MeV and 65 MeV data. The pd data at 13 
MeV in Ref.l2ip are more than 30% smaller than the nd data. A theoretical study of 
the Coulomb effect for the symmetric space star configuration in Ref. 8) shows that 
it is generally very small irrespective of the energy (see Fig. 6 of Ref. 8)). Our results 
at £Ji a b = 13 MeV, 6\ = 02 = 50.5° and (ftu = 120° are located just between the lower 
pd data and the higher nd data, which is very similar to other predictions by the 
meson-exchange potentials. In the other geometrical configurations at 13 MeV, the 
cross sections in the space star 2 case {6\ = 25°, 62 = 50.5°, 4>\2 = 120°) are about 
half of the experimental values and those in the space star 3 case (#1 = 39°, 62 = 
50.5°, 0i2 = 120°) are almost 30% smaller than the experiment. The same situation 
happens in the Faddeev calculations by the Paris potential.'^ (See Figs. 29 and 30 
of Ref. l27p .) Note that we need enough partial waves for the convergence of the 
symmetric space star configurations in particular, which was already pointed out in 
Ref. I]). 

3.6. Coplanar star configurations 

Let us move to the coplanar star configurations in Figs. [7] and [8l The agreement 
between our calculation and the data is satisfactory in general, but some deviations 
still exist. For example, in the first panel of 13 MeV, the new nd datgP are much 
closer to the calculation than the old datap3 but some underestimation still exists 
in the calculation. The underestimation of the cross sections at the minimum point 
5 = 9 MeV of CST2, and also at the np final state interaction peaks around S = 
10 MeV in CST3 and in CST4 are the common feature with the meson-exchange 
potentials. See Figs. 11, 12 and 13 of Ref. l27|) . For 16 MeV reactions, we have given a 
comparison not only for the coplanar star configuration, but also for the intermediate 
star (1ST) configuration, both of which are very similar to the predictions by other 
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models given in Ref. l28|) . For = 22.7 MeV data, the curves are not plotted as 
a function of S but by E2 for the second particle. For this and = 65 MeV 
reactions, we find a large contribution of higher partial waves up to I max = 4. In 
the symmetric backward plane star configuration of 65 MeV, denoted by CST2, the 
original experimental data are shifted to the larger side of S by 3.5 MeV, since the 
starting position of S = does not seem to be the same between our calculation and 
the experiment. 

3.7. Collinear configurations 

The comparison for the collinear configurations are displayed in Figs. M andllUi 
For these configurations, the comparison with the experiment is generally good. The 
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Coulomb force has an appreciable effect to increase the breakup cross sections at the 
collinear point, especially on the low-energy side,® which is much more important 
than the 3N force effect. In the first panel with Ej a t, = 10.5 MeV (COLLI), we find 
some kinematical mismatch of the final state interaction peak at S ~ 10 MeV. For 
COLL2 - COLL5 with E\ & \> = 13 MeV, the small breakup cross sections around the 
collinear points (minimum points) move to better direction to fit the experimental 
data by the expected Coulomb effect. This would also be true for the minimum 
point for E^ = 19 MeV. On the other hand, the breakup cross sections in COLLI 
- COLL4 for E^ = 65 MeV seem to be slightly overestimated. 

3.8. Non-standard configurations 

The comparison with the experimental data for the non-standard configurations 
is shown in Figs. 1111 and 1121 Here we find big deviation from the experimental data 
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in some cases, again a common feature with the meson-exchange predictions. These 
are the NS2 nd scattering of 13 MeV, and the pd scattering of 22.7 MeV and 65 MeV. 
A large number of figures, NS1 - NS9, for 13 MeV are very similar to the predictions 
by the Paris potential in Ref.,27). The huge final state interaction peaks in NS4, 
NS7 and NS8 are very similar to the results by the Malfliet-Tjon potential. In 22.7 
MeV and 65 MeV cases, the calculated results are completely off the experimental 
data. 

§4. Summary 

In this and previous papersf^'^'^ we have applied the quark-model (QM) 
baryon-baryon interaction fss2 to the neutron-deuteron (nd) scattering in the Fad- 
deev formalism for composite particles. The main motivation is to investigate the 
nonlocal effect of the short-range NN interaction in a realistic model, reproducing 
all the two-nucleon properties and yet based on the naive three-quark structure of 
the nucleons. The calculations are carried out by the 15-point Gaussian nonlocal 
potential constructed from fss2, which is accurate enough to reproduce the con- 
verged triton binding energy of fss2 with the accuracy of 15 keV and the NN phase 
shift parameters with the difference of less than o.l .^'^ The potential keeps all 
the nonlocal effects of the original fss2, including the energy-dependent term of the 
QM resonating-group method (RGM). This energy dependence is eliminated by the 
standard off-shell transformation utilizing the square root of the normalization ker- 
nel for two three-quark clusters. It is extremely important to deal with this energy 
dependence properly, since an extra nonlocal kernel from this procedure is crucial to 
reproduce all the elastic scattering observables below E n < 65 MeVp^'^ 

In this paper, we have studied the neutron-induced deuteron breakup differen- 
tial cross sections for the incident energies E n < 65 MeV, and compared them with 
available experimental data and the predictions by meson-exchange potentials. We 
have found that our calculations reproduce almost all the results for the breakup 
differential cross sections predicted by the meson-exchange potentials, including the 
disagreement with the experiment. This feature is probably related with the struc- 
ture of the direct breakup amplitudes in Eqs. ()2- 18j) and f)2- 19j) . First, they are 
constrained by the elastic scattering amplitudes, f(nis>)(tSc)^ m ^ e i m *ial stage. In 
the final stage of reactions, only the half-off shell two-nucleon t- matrix appears owing 
to the energy conservation for outgoing nucleons. The effect of the completely off- 
shell t-matrix therefore appears only at the stage of solving the basic AGS equations 

~(£'S' )J 

for QL 7 , for which the present investigations imply that the difference between 
our QM NN interaction and the meson-exchange potentials is rather minor. On the 
whole, the agreement with the experimental data is fair, but there exist some discrep- 
ancies in certain particular kinematical configurations, which are commonly found 
both for our predictions and for meson-exchange predictions. In particular, the space 
star anomaly of 13 MeV nd scattering is not improved even in our model. There are 
severe disagreement of breakup differential cross sections in some of the non-standard 
configurations. In our model, some overestimations of cross sections are found at the 
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Fig. 12. The same as Fig. 1111 but for other geometrical configurations and energies. The experi- 
mental data are taken from Refs.|27) (nd) for 13 MeV, [23} (pd) for 22.7 MeV, and© (pd) for 
65 MeV. 
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energy E n = 65 MeV. Since these large disagreements can be resolved neither with 
the Coulomb effect nor by the introduction of the 3iV force, systematic studies from 
more basic viewpoints for the AW interaction are still needed both experimentally 
and theoretically. 

In spite of the apparent disagreement between the theory and the experiment in 
some of the breakup differential cross sections, our QM baryon-baryon interaction 
fss2 was very successful to reproduce almost all other experimental data of the three- 
nucleon system without reinforcing it with the three-body force. These include: 1) 
a nearly correct binding energy of the tritonf^ 2) reproduction of the doublet and 
quartet 5-wave scattering lengths, 2 a and 4 apS> 3) not too small differential cross 
sections of the nd elastic scattering up to E n ~ 65 MeV at the diffraction minimum 
points,— 1 ' 4) the improved maximum height of the nucleon analyzing power A y {9) in 
the low-energy region E n < 25 MeV, although not sufficient,^ and 5) the breakup 
differential cross sections with many kinematical configurations, discussed in this 
paper. Many of these improvements are related to the sufficiently attractive nd 
interaction in the 2 S channel, in which the strong distortion effect of the deuteron is 
very sensitive to the treatment of the short-range repulsion of the NN interaction. 
In our QM NN interaction, this part is described by the quark exchange kernel of 
the color-magnetic term of the quark-quark interaction. In the strangeness sector 
involving the AN and UN interactions, the effect of the Pauli repulsion on the 
quark level appears in some baryon channels. It is therefore interesting to study 
Z^-deuteron scattering in the present framework, to find the repulsive effect directly 
related to the quark degree of freedom. Such an experiment is planned at the J-PARC 
facility.^ 
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Appendix A 

Three-nucleon breakup kinematics 

In this appendix, we will discuss the 3A^ breakup kinematics used in this paper. 
We choose the standard set of Jacobi coordinate in the momentum space as a = 3 
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and set 

p=-{kx-k 2 ), q = ^(2k 3 - fei - k 2 ) , K = k 1 + k 2 + k 3 , (A-l) 

with k a being the momentum coordinate of the particle a in the lab system. We 
also choose the z-axis as the direction of the incident particle and assume that the 
magnitude of the incident momentum is qo in the cm system. This implies that 
the incident momentum is always q cm = q$e z and the incident energy is E cra = 
(3H 2 /4M)qo 2 , either the nucleon or the deuteron is the incident particle. Here, e z 
is the unit vector of the z-axis. In the lab system, the incident momentum k^ and 
the energy £j a b are given by 

for nucleon-incident , 

- h h 2 — 2 _ 3 p 

- 2A/ lab ~~ 8M^° ~ 2 cm 




„ „ for deuteron-incident . (A-2) 

- 4j\,/ ft 'lab — 4M^0 — J ™ 

In the following, all quantities in the initial state are expressed by fci a b = |&iab| and 
E cm , which are determined solely by q$. 

In the experimental setup to detect two particles 1 and 2, the three-particle 
breakup configurations are uniquely specified with two polar angles 9\, 2 , and a 
difference of azimuthal angles <f>± 2 = <j>i — 4>2, in addition to the energy 5 discussed 
below. These angles are defined through cos a = (k a ■ e z ) and tan^Q, = (k ay /k ax ). 
We will choose the x-axis such that <p\ = 7r. 16) If k\ = \ki\ and k 2 = \k 2 \ are 
determined from S, k 3 = \k 3 \ and all other angles in the lab system are calculated 
from 

h = [(fclab _ ^1 cos #1 _ ^2 cos 2 ) 2 + {h sin 6*i) 2 + (k 2 sin 9 2 ) 2 



+2k\k 2 sin 9\ sin 62 cos </>i2 



1/2 



COS 6*3 = — (fclab - k\ cos 9\ - k 2 cos 62) , 
h 

sin 9 3 — —— [(ki sin 9\) 2 + (k 2 sin 9 2 ) 2 + 2k\k2 sin 9\ sin O2 cos </>i 2 ] ^ 2 , 
^3 

cos 63 = - — : — ; - (ki sin 9\ + k 2 sin 6 2 cos 4>\ 2 ) , 
k 3 sin 6I3 

sin <j) 3 = — - — : — — /c 2 sin 2 sin 12 , (A-3) 
h sm 6I3 

with fa = ir and 02 = 7T — <Pi2- Once all the momentum vectors in the lab system 
are determined, the momentum vectors in the cm system are easily determined from 
the relationship, p = (l/2)(fci — k 2 ) and q = (2/3)fci a b — k\ — k2 for a = 3, and 
their cyclic permutations of (123). The basic magnitude of q = \q\ is obtained from 
k 3 in Eq. ()A-3p by a simple replacement of fci a b with (2/3)&i a b: 

q = [((2/3)fc Iab - fei cos 6»i - k 2 cos 6> 2 ) 2 + (&i sin 6>i) 2 + (/c 2 sin 9 2 ) 2 

+2kik 2 sin 6>i sin 9 2 cos 0i 2 ] 1/2 . (A-4) 
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The two-nucleon momentum p = \p\ is determined from the energy conservation in 
the cm system: 

E = E cm + s d = ^(p 2 + p 2 ^) , (A-5) 

where e d = -(3h 2 /4M)n d 2 is the deuteron energy. It is convenient to use the 
threshold momentum qu = sj qo 2 — K d 2 for the deuteron breakup, by which we 
find 



p = \1\(<im 2 -q 2 ) =Po ■ (A-6) 



The angles of p and q are obtained from 

cos 6p = — (k\ cos #i — k2 cos #2) , 
2p 

cos = — (k\ sin 9\ cos <f>i — k2 sin O2 cos fo) 

2psm Up 

sin (f>„ = — (k\ sin 6\ sin d>\ — k% sin 62 sin 62) 

2psm dp 

1 (2 

COS d g = — -^lab — k\ COS 6*i — /C2 COS 6*2 

(/ci sin 9\ cos </>i + k2 sin 6*2 cos (^2) 



q sin „ g 

sin a = — (k\ sin 6\ sin <^>i + /C2 sin 62 sin ^2) . (A-7) 

gsin # g 

In order to determine k\ and /C2 from S, we start from the energy conservation 
in the lab system: 

h 2 

Ex^ + e d = — (k l 2 + k 2 2 + kz 2 ) . (A-8) 

We rewrite this as 

h 2 + h 2 + ^1^2 cos 0i2 - fci a b(^i cos 6»i + A; 2 cos 6> 2 ) + A = , (A-9) 
where cos #12 = (k% ■ k2), and we have defined 

M 1 o f i K d 2 f nucleon-incident 

A = ~la (^ab + Erf) + TT^lab 2 = S t , 2 , o for S 

" I I^iab + I'm [ deuteron- incident 

(A40) 

We rotate k\-k2 plane by 45° and parametrize the ellipse with an angle 6. In this 
process, it is convenient to express #12 by 6q, which is defined by 



Arccos -cos 6 U . (ATI) 
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Here, Arccos implies the principal value between and it. Note that 9q changes from 
7r/6 to 7r/3 for the change of 9 12 from to n. If we use this 9q, the solution of Eq. 
(|A 9|) is parametrized as 

J fciab_4_ ^_4c(0-0 o ), (A-12) 

2^/2 sin 29 V u; ' 2\/2 sin 20 O V 7 V ' 



with 



and 



C{6) = cos 9 + C , S(0) = sin (9 + 29 Q -vr/2) + 5 , (A-13) 



a/2 cos 0i — cos 29q cos #2 c V% cos 6*2 — cos 29q cos #i 
C ° = T" ^29~ ' So = li ^W • (A ' 14) 

In Eq. (|A-12j) . we have defined 

A _ /of ( cos g i) 2 + ( cos g 2) 2 ~ 2 cos 29 cos fli cos fl 2 4Z\ y 

I (sin2fl ) 2 fc lab 2 J • [ ' j 

We measure the arc length S in the E1-E2 plane counterclockwise, starting from 
a certain starting point 9 s t- The expression S(9) is obtained by integrating dS = 



^/{dE l Y + (dE 2 f: 

- 2 fclab 2 / A N 2 

2M V2 sin 26»o y ./« 

with 



/(0) = y (C(0) sin #) 2 + (5(0) cos (0 + 2# - ^/2)) 2 . (A-17) 

Note that S{9) is a monotonically increasing function of 9, satisfying S(9 st ) = and 
S{9 + 2n) = S{2tt) + S (9). 

For the practical calculation, we first discretize the integral region of Eq. (|AT6p 
into small intervals by 

e it = Bt + iijf , (A- 18) 

with n = 1 - N . The number of discretization points N is typically N = 50. We use 
third order spline interpolation for 9, 

N 

f(9-9 ) = Y,St ) (e)f(9 u -9 ) , (A-19) 

with 

3 

Si K) (0) = E a » (m) ( 9 ~ B «T for 9 G [^-1 A] • (A-20) 

m=0 
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If we integrate Eq. (|A-16|) over 8 from 9 st to 8^ by using Eqs. f)A19|) and (|A-20|) . we 
can carry out the ^-integral analytically and obtain 



h 2 h - 2 



Hab 



A 



2M V 2 sin 28, 



m=0 



£(-i) m+1 ^(f ) EE< (m) /(^-^o) . (A-2i 



U=l K=l 



To obtain the angle 8 from the arc length 5*, we again use the spline interpolation 
technique, 



N 



0(S) = J> M (S)0(5„) 



(A-22) 



Here, 5^(5) is the third order spline function for the mesh points {S^} = [So, Si, S2, 
■ ■ ■ , Sn] with So = 0. From Eq. (|A-22j) with 8(S^) = 8^, we obtain 



8{S) 



2vr 



N 



(A-23) 



where Eq. (|A-18|) and ^2u=i S^(S) = 1 are used. We therefore only need to calculate 
the sum in Eq. (|A-21|) and prepare the coefficients of the spline interpolation for the 
mesh points {S^}. 

The starting angle 8 is selected as follows. Let us first consider the nucleon- 
incident reaction. In this case, it is convenient to define the angle 8^ throughP 



cos 



-) 



2y/A 2 K d 



(A-24) 



hah s/3 qo 

If we assume /c2=0 in Eq. (]A-9p . we find that there are two non- negative solutions 



(±) _ hab 



cos 81 ± J (cos 8i) 2 - (cos 8(-)) 2 



(A-25) 



only when 81 < 8^ \ We choose the larger value k^ as the starting point to measure 
S. In this case, we can easily find that the corresponding 8 value is given by 



8 s t = Arccos Sq — 80 



(A-26) 



where So is given in Eq. (|A"T4j) . When 81 > flH, we have two cases. In the case of 
82 < 8^~\ we choose ki = and 



.(-) _ 



as the starting point with 



cos 8 2 - V (cos 8 2 ) 2 - (cos 8(-)) 2 



(A-27) 



6 s t = Arccos Cq + 80 — it ■ 



(A-28) 
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In the case of 62 > Q^~\ the ellipse does not cross over either k\ or k<i axis. We 

therefore use the smaller value of k\ = k^ as the starting point. This condition 
yields 

COS &2 



. 1 cos 9\ 

Arccos ( — = 

A y/2 sin O 



+ 7T 



(A-29) 



and 



.(-) 



.(-) 



4 cos 



cos 9\ + cos 62 
2 cos 6*o 



'cos 0i + cos 2 \ 2 4Z\ 



2 cos 0o 



(A-30) 



When the deuteron is an incident particle, there is no crossing point across either 
the fci-axis or the ^-axis. We follow the definition of the Correll et al.'s paper 
that discusses the deuteron incident reaction around the collinear configurations, 
choosing the collinear point as the starting point to measure S. The collinear point 
C is defined as the configuration with q = (2/3)fei a b — {k% + £2) = 0. To find the 
corresponding st = C , we assume 0i2 = vr (0i = n, 4>2 = 0) and the x-z plane 
as the reaction plane, just as the experimental setup. Under this assumption, the 
magnitude q 2 is expressed as 



fciab I + (&l sin 0i - k 2 sin 2 ) 



q~ = ( ki cos 0i + k2 cos 02 
resulting in the two conditions 

k\ cos 0i + k2 cos 02 = TT^iab , ^i sin 0i = fc 2 sin 02 • 

If 0i = 02, we assume k\ = A; 2 and find C = tt/2 with 

cos 01 



jj^lab 



fei = ko 



This corresponds to the k^ = k^ 1 case of Eq. (|A-30j) with the opposite sign for the 
second term. In the general case, the three conditions of Eq. (|A-32p and the energy 
conservation in Eq. (|A-9j) are not always simultaneously satisfied. (Note that we 
only need two conditions to determine k\ and fo.) We take the following procedure 
to determine 9 r . Let us use the notation 



fclab 

4 

.(+) 



+ 



A 



(cos O ) 2 v 7 ^ cos O 



(A-31) 



(A-32) 



(A-33) 



COS 01 — COS 02 



V2 



sm t>Q 



(A-34) 



to simplify the expressions. We have a 2 + b 2 
angle a byEE 



cos 0i + cos 02 

\fi COS 0o 

A 2 + (8Z\/fci ab 2 ) and define a new 



b 



cos a 



Va 2 + b 2 



sm a 



Va 2 + b 2 ' 



(A-35) 



*' By using a, Co and So in Eq. (IA-14|) are expresses as Co 
VI + e 2 sin (a + 9 - n/2) with e = (2V23/fciab^)- 



vT + e 2 cos (a — So) and So = 
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We consider, 1(6) = k\ cos 9\ + k 2 cos &2, ^ a function of 9, by using Eq. (|A-12p and 
others. We find 



Y 

1(9) = ki cos 9i + k 2 cos 9 2 = -MabV a 2 + b 2 A cos (9 - a) + \J a 2 + 6 2 

The crossing point with 7(0) = (2/3)/ci a b is found only when the condition 



1 

A 



1 



3 Va 2 + £> 2 

is satisfied. The solution 9 = 9 C is found as 



Va 2 + b 2 



< 1 , 



9 r = a ± Arccos 



1 



1 



\J a 2 + 6 2 



(A-36) 



(A-37) 



(A-38) 



.3 Va 2 + & 2 

and a unique point is determined when the equality is satisfied in Eq. (|A-37p . Next, 
we examine the condition, k\ sin 9\ = k 2 sin 9 2 , is satisfied or not for the two so- 
lutions of Eq. (|A-38p . The one satisfying this condition is the collinear point with 
q = from Eq. (|A-3ip . If both solutions satisfy the condition, we choose the smaller 
one for 9 C . If neither of the solution satisfies the condition, there is no exact collinear 
point. In this case, we minimize Eq. (|A-3ip with respect to 9 in the interval bounded 
by the two solutions of Eq. (|A-38|) . 



Appendix B 

Isospin factors for the breakup amplitudes 



In this appendix, we extend the definition of the spin factors 



S ^ S S z fS z 



(-i) l+ °xl 



for a 



1 

2 
3 



(B-1) 



to the isospin factors and calculate the matrix elements (r T \Or^\r r ) in Eq. ()2- 13j) . 
The most convenient definition of the isospin factors is probably 

(r T \of\r T ) = ((^)|r 2 |(p ( \ 2 3 ) ) 3 - a o r ( J Pf 123 )) /3 l(^)^ 2 ) 



2/2" 

\i+t Yj( a ^) 
t,t 

(_l)l+tvT(af3) 
^ ' t,t 



( (-l) 1+t Xl { 



for f3 — a = < 



1 

2 in (mod 3) 
3 



(B-2) 



which yields the results in Eq. (|2-14|) . If we set O t = 1, all the factors X~^ a ^ in 



Eq. (|B-2p are reduced to X~ since Of = (P[ 123) f- a . Here, X 1 / 2 i s the common 



t.t 



(123) ' 



matrix with the spin factors given by^ 



X ls< 



1 


V3 


2 


2 




1 


2 


2 



3 

X 2 . 




1 



(B-3) 
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In Eq. ()B-3p . the upper row (the left-most column) corresponds to s = (s' = 0) 
and the second row (the right-most column) corresponds to s = 1 (s' = 1). 

In order to calculate X~^\ we only need (t\0 T \t) with \t) = \{t\)\T z ), since 
P(i23) does not change the total isospin T = 1/2. Furthermore, O r in Eq. (|2T0p 
are expressed by the symmetric isospin operator t z = (t z (1) + r 2 (2))/2 and the 
antisymmetric operator = ( T ^(l) — T ^(2))/2 for the two-particle states 7^(1,2). 
The former does not change t = or 1, while the latter flips the isospin value. For 
pp or nn, the non-zero matrix elements are only for t = t = 1, but np and pn 
contains t z . However, we only need to calculate the sum of np and pn contributions. 
Furthermore, Eq. (|2-5p tells us that np and pn give the same contribution owing to 
the permutation operator (1 + P). We can therefore assume O t in Eq. (|F3-2[) as 

OPP = h z (l + t z ) , O nn = U z {-l + t z ) , O pn + O np = 1 - t 2 z . (B-4) 

If we decompose t z into the rank 0, 1, and 2 tensors as t 2 z = (l/3)t 2 + ^2/?>[t i]g , we 
immediately find that the rank 2 tensor does not contribute since the total isospin in 
our case is T = 1/2. We can therefore replace t 2 in Eq. (|B-4p with (l/3)i 2 , resulting 
in 

{t\O p P\t) = (t\o nn \t) = s lt s t ^ , 

(tpnt) + (t\O np \t) = 5 lt (Vo + ^,i) . (B-5) 

It is convenient to introduce the isospin projection operators Pq = (1 — r(l) • r(2)) /4 
and Pi = (3 + r(l) • r(2))/4, and define 



(123)/ - 1 T V' t (123)' 
-l) 1+ *(i|(^(l 23 )) 3 ^(^23)) /3 | / ) 4 

\* 1(^(123)) "-^(^(m)) 01 !* 



' 1 



2 in (mod 3) , 
3 

(B-6) 



for t = and 1. Since P T is expressed as P T = \t){t\ in our model space, the matrix 
elements {t\{P( m ?- a Pr{P{ m f\t) = ft(P{ 123) ) 3 - a \r) (r\(Pf m f\t) can be easily 
calculated from Eq. (|B-ip . The correspondence 

qpp, O nn ~ -Pi , O pn + O np ~ P + -Pi , (B-7) 
3 3 

from Eq. (|B-5p yields 

2 1 

jfpp^ X nn = -X 1 ^^ , X p ™ + X np = X°( a ^ + -X 1( - a ^ , (B-8) 

3 3 

in the matrix form. Here X T ( a/3 ^ with r = and 1 are given by 
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A -(H22) _ I 4 \ X Q(33) / <■< 



V3 3 / ' V o o 

4 4 



X 0(12) = X 0(21) 



_ 1 

4 4 

4 4 



(r = 1 factors) 

A -UJU I 3 -T | x i (22 ) = / i ^ I X l(33) ^00 



va i I ' \ %/i i / ' V 01 

4 4/ \ 4 4 



X l(12) = X l(21) 



3 _VJ 

4 4 

yj i 

4 4 







-f 
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